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Abstract. New SiO multi-transition millimetre line observations of a sample of carbon stars, including J = 8— »7 observations 
with the APEX telescope, are used to probe the role of non-equilibrium chemistry and the influence of grains in circumstellar 
envelopes of carbon stars. A detailed radiative transfer modelling, including the effect of dust emission in the excitation analysis, 
of the observed SiO line emission is performed. A combination of low- and high-energy lines are important in constraining the 
abundance distribution. It is found that the fractional abundance of SiO in these C-rich environments can be several orders of 
magnitude higher than predicted by equilibrium stellar atmosphere chemistry. In fact, the SiO abundance distribution of carbon 
stars closely mimic that of M-type (0-rich) AGB stars. A possible explanation for this behaviour is a shock-induced chemistry, 
but also the influence of dust grains, both as a source for depletion as well as production of SiO, needs to be further investigated. 
As observed for M-type AGB stars, a clear trend that the SiO fractional abundance decreases as the mass-loss rate of the star 
increases is found for the carbon stars. This indicates that SiO is accreted onto dust grains in the circumstellar envelopes. 
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1. Introduction 

In addition to dictating the evolutionary time scale, stellar 

winds from late-type stars are thought to be important for the 

J '• II 1 

enrichment of heavy elements in the ISM JGustafss on 2004). 

The stellar wind carries the results of internal nuclear pro- 
cesses, activated during evolution along the asymptotic giant 
branch (AGB), and thus contributes to the chemical evolution 
of galaxies. 

Red giant stars located on the AGB divide into two chemi- 
cally distinct groups; carbon stars (C/0> 1), which show a more 
rich chemistry than the M-type (0-rich) stars (C/0<1). The 
molecular- and grain-type setups in the circumstellar envelopes 
(CSEs) around these mass losing stars are to a large extent de- 
termined by the C/O-ratio of the central star. Of paramount im- 
portance is to know how much of the elements that are eventu- 
ally locked in dust particles and how much remains in the gas 
phase. 

Observations of SiO line emission is particularly inter- 
esting because it is potentially a useful probe of the forma- 
tion and evolution of dust grains in CSEs, as well as CSE 
dynamics: hinted a t from early interferom etric observations 
llLucasetaI.lll992t ISahai&Bieging|ll993l) . A major survey 
of SiO emission from M-type AGB stars was performed by 



iGonzalez Delgado et alJ ( l2003h . Supplemented by a detailed 
numerical modelling of the detected radio Une emission, they 
found that the circumstellar SiO fractional abundance (relative 
to H2) for irregular and semi-regular variables in their sam- 
ple had a median value of 5x10"^, a factor of about ten lower 

I . II I 
than expected from theory (see, e.g., review by Millar 2003|). 

For the high-mass-loss-rate Miras in their sample Gonzalez 
Delgado et al. found a distinct division into a low-abundance 
group (on average 4x10"^) and a high-abundance group (on 
average 7x10"^). They concluded that this division is not an 
effect of the modelling, but they found no obvious astrophysi- 
cal explanation. 
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Interestingly, Gonzalez Delgado et al. (2003) found a trend 
of decreasing circumstellar SiO fractional abundance with in- 
creasing mass-loss rate of the star, interpreted as an effect of 
increased adsorption of SiO onto dust grains with increasing 
mass-loss rate, i.e., the density of the envelope. This claim 
has been further corroborated by the recent interferometric SiO 
observations of the two low-mass-loss-rate M-type AGB stars 
R Dor and L^ Pup performed bv Schoier et al. (2004) . The exci- 
tation analysis of these SiO (v=0, 7=2— >1) interferometric ob- 
servations requires the fractional abundance of SiO to be high, 
a;4x 10"^, within w 1 x 10'^ cm as suggested by stellar atmo- 
sphere models. Beyond this radius the SiO fractional abun- 
dance drops by about an order of magnitude and at a; 3 x 10'^ 
cm photodissociation destroys the SiO molecules. Thus, there 
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Table 1. Integrated (lobs - j T^h dv) line intensities in K km s"' for the new observations of SiO (v = 0, 7 — > 7 - 1) line emission. 



Source" 



OSO 



NRAO 



SEST 



JCMT APEX 



1 3->25^4 2-^13->25-^4 



LP And 

RVAqr 

UCam 

SCep 

VCyg 

RFor 

V821 Her 

CWLeo 

RLep 

RWLMi 

V384 Per 

WPic 

RVol 

IRAS 07454-7112 

IRAS 15082-4808 

IRAS 15194-5115 




20 



" A colon (:) marks a low S/N detection. ^^^^_^^^^ ^_^ 

* Improved S/N-ratio compared to previous detection bv lOlofsson et alj il998l) . 



3.87 



2.62 
4.97 

2.73 



0.94: 

5.12 
8.51 
20.4 



exists strong indications that the circumstellar SiO line emis- 
sion carries information on the region where the mass loss is 
initiated, and where dust formation takes place. 

In this paper we present new (sub-)millimetre line observa- 
tions of SiO for a sample of carbon stars. The observations are 
supplemented by a detailed non-LTE radiative transfer analysis 
in order to obtain reliable circumstellar molecular abundances. 
The results are then compared to predictions from available 
chemical models and to the abundance estimates of SiO in M- 
type (oxygen-rich) AGB stars. 



2. Observations 

The observations were performed during 1994-1998 using the 
Onsala 20 m telescope' (OSO), the Swedish-ESO submillime- 
tre telescope^ (SEST), the NRAO 12 m telescope at Kitt Peak, 
in October 2003 using the JCMT telescope^, and in September 
2005 using the APEX telescope"^. The observed spectra are pre- 
sented in Figs[2&|2and velocity-integrated intensities are re- 



Table 2. Telescope data. 



' The Onsala 20 m telescope is operated by the Swedish National 
Facility for Radio Astronomy, Onsala Space observatory at Chalmers 
University of technology 

^ The SEST was located on La Silla, Chile and operated jointly by 
the Swedish National Facility for Radio Astronomy and the European 
Southern Observatory (ESO). 

^ Based on observations obtained with the lames Clerk Maxwell 
Telescope, which is operated by the Joint Astronomy Centre in 
Hilo, Hawaii on behalf of the parent organizations PPARC in the 
United Kingdom, the National Research Council of Canada and The 
Netherlands Organization for Scientific Research 

"* This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fr Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory 



Transition 


Frequency 
[GHz] 


Telescope 


'/mb 


["] 


7 = 2^1 


86.847 


SEST 


0.75 


57 






OSO 


0.65 


44 


7 = 3^2 


130.269 


NRAO 


0.80 


49 






SEST 


0.68 


38 


7 = 5-^4 


217.105 


NRAO 


0.55 


29 






SEST 


0.50 


24 






JCMT 


0.70 


23 


7 = 8^7 


347.331 


APEX 


0.70 


18 






JCMT 


0.62 


14 



ported in Table [2 The intensity scales are given in main-beam 
brightness temperature scale (Tmb)- 

The SEST, OSO and JCMT observations were made in a 
dual beamswitch mode, where the source is alternately placed 
in the signal and the reference beam, using a beam throw of 
about 11' (SEST and OSO) or 2' (JCMT). This method pro- 
duces very flat baselines. At the APEX 12 m telescope the ob- 
servations were carried out using a position-switching mode, 
with the reference position located +2' in azimuth. The raw 
spectra are stored in the T'^ scale and converted to main-beam 
brightness temperature using r^b = T^/ri^b. T* is the antenna 
temperature corrected for atmospheric attenuation using the 
chopper-wheel method, and /^mb is the main-beam efficiency. 
Regular pointing checks were made on SiO masers (SEST and 
OSO) and strong CO sources (APEX) and typically found to 
be consistent within x; 3". 

At the NRAO 12 m telescope the observations were carried 
out using a position-switching mode, with the reference posi- 
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Fig. 1. New observations of SiO line emission. 
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Fig. 2. New observations of SiO 7 = 8— > 7 line emission using the APEX telescope. The velocity resolution is l.Okms ' in all 
the spectra, except for WPic which has been reduced to a resolution of4.0kms'. 



tion located +10' in azimuth. This is the preferred observation 
mode for spectral line observations at the 12 m telescope. The 
raw spectra, which are stored in the T^ scale, were converted 
using T^b-Ty?]*^, where rj*^, is the corrected main-beam effi- 
ciency. The T^ scale is related to T'^ through T^-Tyrjfss, where 
;7fs,s is the forward-scattering and spillover efficiency. Regular 
pointing checks were made on strong continuum sources and 
typically found to be consistent within x 5". 

The adopted beam efficiencies, together with the FWHM 
of the main beam (0mb), for all telescopes and frequencies are 
given in Table|2] The uncertainty in the absolute intensity scale 
is estimated to be about ±20%. 

The data was reduced in a standard way, by removing a low 
order baseline and then binned in order to improve the signal- 
to-noise ratio, using XS^. 

In addition to the new data presented here we have also 
used Si O Une intensities re ported by Olofsson et al. (1982 , 
1998) & IWoodsetal1 l l200j^hr./=2^n .'Ruiarra baletalJlll994h 
(7=2^1 and 3^2), and iBie g in get al.. (.20001) (7=5^4 and 
8— >7). In total, 19 carbon stars have been detected in SiO line 
emission. 

These stars detected in SiO millimetre line emission pro- 
vide our sample and are listed in Table|3l 

3. Molecular excitation analysis 

3.1. Radiative transfer model 

The CSEs are assumed to be spherically symmetric, produced 
by a constant mass loss rate (M), and to expand at a constant 
velocity (ve). In order to determine the molecular excitation 
in the CSEs a detailed non-LTE radiative transfer code, based 



^ XS is a package developed by P. Bergman to reduce and analyse 
a large number of single-dish spectra. It is publically available from 
[ftp : //yggdrasil . oso . Chalmers . se . 



on the Mon te Carlo method, wa s used. The code is described 
in detail in ISchoier & Olofssonl ll200ll) and has been bench- 
marked, to high accuracy, agai nst a wide var iety of molecular- 
line radiative transfer codes in Ivan Zadelhoff et aL C200Z) and 
van der Tak et al. (in prep.). 

ISchoier & Olofssonl ( 1200 la set out to determine the circum- 
stellar physical properties, such as the density, temperature, and 
kinematic structures, of a large sample of carbon stars based 
on radiative transfer modelling of millimetre CO line obser- 
vations. Thes e models, and subseque n t mod els presented in 
ISchoier et all (§002) and IWoods etaD (l2003l) . form the basis 
of the underlying physical structure and are used as input to 
the SiO excitation analysis. In a similar way new circumstel- 
lar models where obtained for V821 Her (IRAS 18397H-1738) 
and VI 965 Cyg (IRAS 1 9 32 1-H2757) usin g CO milli metre line 
observations published in lMargulis et al.l ( ll990i) . iNyman et alJ 
|l992), and from the JCMT data archive^. 

The local line width is assumed to be described by a 
Gaussian and is made up of a micro-turbulent component and a 
thermal component. The micro-turbulenc e is assumed to have 
a width of 1 .0 km s ' dSchoier et al.l2004l) . This is significantly 
higher than the value of 0.5 km s' adopted in the CO mod- 
elling performed by Schoier & Olofsson (2001). However, as 
pointed out by Schoier & Olofsson the exact value adopted 
for the micro-turbulence only has a minor effect on the de- 
rived mass-loss rates, as long as it is much lower than the ex- 
pansion velocity of the wind. The (sub)millimetre SiO obser- 
vations typically probe regions of gas closer to the star than 
does the CO sub (millimetre) observat ions. Observations of IR- 
absorption lines dMonnier et al.l2000l) indicate that most likely 
there is a gradient in turbulent velocity, with regions closer to 
the photosphere being more turbulent. Thus, it is not unreason- 
able to adopt a higher turbulent velocity when modelling the 



|http : //www ■ j ach . hawai i . edu/3 CUT/ archive/ 1 



Table 3. Model results. 
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" U Cam has a detached shell that complicates the analysis of the present-day mass-loss characteristics (for details see lSchoier et al 12005 ah . 
A colon (:) marks an uncertain abundance estimate. 



SiO line emission. The thermal line broadening is calculated 
from the derived kinetic temperature structure. 

Since the publication of ISchoier & Olofssonl ( 1200 ll) new 
colhsional rat e coefficients of CO have been calculated by 
iFlowen (1200 ll) . These new rate coefficients were adopted here 
and extended to include more energ y levels as well as extr apo- 
lated in temperature, as described in lSchoier et alJ (l2005bl) . An 
ortho-to-para ratio of 3.0 was adopted when weighting together 
collisional rate coefficients for CO in collisions with ortho-H2 
and para-H2. In addition, dust emission has been included in 
the excitation analysis for CO in the same fashion as for SiO 
described below. Moreover, new CO observations of / = 3 — > 2 
line emission using the JCMT and APEX telescopes (Schoier 
& Olofsson, in prep.) has been included in the modelling. In 
all, these new modifications to the excitation analysis of CO 
millimetre line emission have not alt ered the mass-loss rates 
derived bv*S choier & Olofssonl (120011) by more than 20%, i.e., 
within the calibration uncertainty of the observations. The new 
mass-loss rates are reported in Table|3 

The excitation analysis includes radiative excitation 
through the first vibrationally excited (v=l) state for both 
CO at 4.6 ;um an d SiO at 8|/m. Relev ant molecular data 
are summarized in ISchoier et alJ (l200 .5b) and are made pub- 
licly available through the Leiden Atomic and Molecular 
Database (LAMDA)^. The majority of sources in this study 
are intermediate- to high-mass-loss-rate objects where thermal 



dust emission provides the main source of infrared photons 
which excite the v = 1 state. 

The addition of a dust component in the Monte Carlo 
scheme is straightforward as described in lSchoier et alJ ( l2002l) . 
At the wavelengths of importance here (/l>4.8jum) scattering 
can be neglected and only emission and absorption by the dust 
particles need to be considered. The number of model photons 
emitted per second by the dust within a small volume of the en- 
velope. Ay , and over a small frequency interval, Av, can then 
be written as 

mr) = ^K,pA(r)BATi{r)]^V(r)^v, (1) 

nv 

where Ky is the dust opacity per unit mass, pd the dust mass 
density, and By the Planck function. T^ir) is the temperature 
of the thermally emitting dust grains. Typically, the frequency 
passband corresponds to 3xve and it is centered on the line rest 
frequency vq. 

The model photons emitted by the dust are released to- 
gether with the other model photons. The additional opacity 
provided by the dust. 



(2) 



' .http : //www ■ strw . leidenuniv . nl/~moldata| 



Ty-Ky I piir) dr, 

is added to the line optical depth. 

3.2. SED modelling 

The dust-temperature structure and dust-density profile are ob- 
tained from detailed radiative transfer modelling using Dusty 
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dlvezic & Elitzun 19971) . Here we are only interested in the po- 
tential of dust grains to affect the excitation of molecules. A full 
treatment of dust radiative transfer coupled with a dynamical 
model for the wind, in order to independently measure physi- 
cal properties of the CSE such as, e.g., the mass-loss rate, is 
beyond the scope of this paper and left to a future publica- 
tion (Schoier & Olofsson, in prep.). In the dust modelling the 
same basic assumptions are made as for the gas modelling, i.e., 
a spherically symmetric envelope expanding at a constant ve- 
locity. Prompt dust formation is assumed at the inner radius, 
ri, of the CSE. An iorphous carbon dust grains with the optical 
constants given in Suh (2000) are adopted. For simplicity, the 
dust grains are assumed to be of the same size (a radius, ad, 
of 0.1 //m), and the same mass density (ps = 2.0gcm"^). The 
corresponding dust opacities, Ky, were then calculated from the 
optical constants and the individual grain properties using stan- 
dard Mie theory (Bohren & Huffman 1983). 

In the modelling, where the SED provides the observa- 
tional constraint, the dust optical depth specified at lOyum, 
Tio, and the dust condensation temperature, rd(ri), are the 
adjustable parameters in the ;t'^-analysis. The model SED 
only weakly depends on the other input parameters which 
are fixed at reasonable values. The SED typically con- 
sists of JHKLM photometric data (' KerschbaumI 1 1 9991 and 
Kerschbaum priv. com.), IRAS fluxes, and in some cases sub- 
millimetre data (Groenewegenet al. 1993). The total luminos- 
ity of the source is obtained from the pe riod-luminosity rela- 
tion o flGroenewegen & Whitelockl(ll996l) . The distance is then 
obtained from the SED fitting. The derived parameters are re- 
ported in Tableland examples of best-fit models are shown in 
Fig.|3]for the new sample sources V821 Her and V1965 Cyg. 
A more detailed description of the procedure as well as other 
examp les for the sample stars can be found in ISchoier et alJ 
1I2OO2I) . 

3.3. SiO envelope sizes 

In order to determine accurate abundances it is im- 
portant to know the spati al extent of the molecules. 
lOonzalez Delgado et al J (120031) showed that it is possible to put 
constraints also on the size of the SiO line emitting envelope 
from single-dish observations alone, given that enough transi- 
tions are observed (in practice four or more). In the present 
sample six carbon stars (CW Leo, R For, R Lep, RW LMi, 
S Cep, and V Cyg) have been observed in four or more tran- 
sitions. For these sources the size of the SiO emitting region 
(re) and the inner fractional abundance of SiO (/o) are varied 
simultaneously. The abundance distribution is assumed to be 
described by a Gaussian 



/(r) = /() exp|-(- 



2\ 



(3) 



where / = «(SiO)/«(H2), i.e., the ratio of the number density 
of SiO molecules to that of H2 molecules. In circumstellar en- 
evelopes such as these H is expected to be mainly in molecular 
form. 
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Fig. 3. Best-fit models, derived from the dust radiative transfer 
modelling (parameters given in Table|5}, of the observed SEDs 
for V821 Her and VI 965 Cyg . 

The best fit model is found by minimizing the total x^ de- 
fined as 



/^'to 



/=1 



v^mod ^obs7 



(4) 



where / is the integrated line intensity and cr the uncertainty 
in the measured value (usually dominated by the calibration 
uncertainty of ±20%), and the summation is done over A^ in- 
dependent observations. The sensitivity of the line emission to 
variations of the two adjustable parameters r^ and /o is illus- 
trated in Fig. 0] As can be seen reasonable estimates can be 
found for all four objects, also the quality of the best-fit models 
are good, typically ;\f^g^ -X^t/i^ - 2) ~ 1. The Icr range of ac- 
ceptable values for r^. and /o are reported in Table 0] However, 
note that the range of acceptable values does not follow a rect- 
angular shape as seen in Fig.0] 

Given the relatively low number of carbon stars where r^ 
can be determined we have assumed it to scale with density as 



M 
log re = 19.2 H- 0.48 log I — 



(5) 
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Fig. 4. x^ maps showing the sensitivity of the excitation analysis to the adopted SiO fractional abundance (/o) and envelope size 
(re). Contours are drawn at the 1, 2, and 3cr levels. Also indicated is the lowest reduced;^^^ (x^^^ - xloili^ ~ 2)) in the map. 



Table 4. SiO abundance distribution. 



JO 



Source 



[cm] 



N xL 



SCep 


1.4 ±0.8x10 


VCyg 


1.3 ± 0.7 X 10 


RFor 


2.2 ± 1.8x10 


CWLeo 


2.8 ± 1.1x10 


RLep 


3.1 ± 1.4x10 


RWLMi 


6.0 ± 4.0 X 10 



5 4.4 ±1.6x10'^ 5 0.6 

5 3.6±1.3xl0'5 6 2.0 

5 5.4 + 3.6x10'^ 4 0.3 

'-' 1.9 + 0.7x10"^ 9 0.8 

'-« 4.4±1.6xl0'' 4 0.9 

7.7 + 2.5x10'^ 5 1.1 



" The abundance distribution is assumed to be Gaussian (see Eq.|3j 

where M is the mass-loss rate and Vg is the expansion 
velocity of the wind. This scaling law was found by 
lOonzalez Delgado et alJ 1120031) for a larger sample of M-type 
AGB stars. As shown in Fig.|5lthis scaling law can, within the 
uncertainties, account for the estimated envelope sizes derived 
for the carbon stars in our study. 

The radial distribution of SiO molecules in the gas- 
phase is most probably determined by a combination of con- 
densation and photodissociation processes. Thus, one can 
imagine an initial (pre-condensation) SiO abundance deter- 
mined by the stellar atmosphere chemistry, possibly modi- 
fied by non-LTE chemical processes due to shocks in the 
inner part of the wind ou t to a few stellar radii (e.g., 
IWillacv & Ch erchneff 1998). In terms of classical evapora- 
tion theory (see G onzalez Delgado et al. 2003. and references 
therein) the region over which the SiO abundance is expected 
to change significantly due to condensation onto dust grains 
is confined within approximately 10'^ cm. For intermediate- to 




Fig. 5. SiO envelope sizes (re) estimated from the excitation 
analysis plotted against a density measure (M/ve). The solid 
line is the relation found for a larger sample of M-type AGB 
stars (Eq.|5|i. 



high-mass-loss-rate objects this region is significantly smaller 
than the post-condensation region (see Fig. |5}, which is lim- 
ited by photodissociation, and it does not contribute notably 
to the observed single-dish fluxes. Thus, a single component 
is adequate for describing the SiO abundance structure in this 
case. It should be pointed out that this also means that no con- 
straints can be put on the pre-condensation abundance from the 
present observations. This requires high spatial resolution ob- 
servations t.Schoier et al. 2004) . In the case of low-mass-loss- 
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rate objects the post-condensation region is smaller, but, since 
the freeze out is less effective, the contrast between pre- and 
post-condensation SiO abundances is small. Hence, a single 
component description is adequate also in this case. 

3.4. SiO abundances 

Using the scaling relation in Eq. |5] it is possible to estimate 
the characteristic size of the envelope (re) containing SiO 
molecules for each of the sample sources. Once the size is 
known only one free parameter remains, the fractional abun- 
dance of SiO (/o). The derived sizes and abundances are re- 
ported in Table 13 The fits are generally excellent with reduced 
X^ values of a!l-2. As examples we show the best-fit models 
for R For and R Lep in Fig. |6l 

From the modelling it is noted that the SiO J = 3 —> 2 line 
observ ed using the IRAM 30 m telescope by .Buiarrabal et alJ 
( 1 19941) is consistently lower by about a factor of two compared 
to model predictions. A possible explanation for this is pointing 
offsets on the order of x 8". Due to the larger beam the J = 
2 — > 1 observations are less affected, typically on the same 
level as the absolute calibration uncertainty ^ 20%. However 
a calibration error can not be ruled out. In the present analysis 
the SiO 7 = 3 — > 2 line intensities reported by Buiarrabal et alJ 
ll 19941) were not included, except for UU Aur which is only 
detected in this line. 

Making an error estimate of the derived abunda nce is dif- 
ficult as discussed in lOonzalez Delgado et al.l (l2003h . We have 
also tested the effect of the assumed abundance distribution by 
replacing the adopted Gaussian distribution (Eq.|3 by an expo- 
nential decline of the abundance, f(r) - fi, exp(-r/re), as would 
be the case for a simple photodissociation model without any 
dust shielding. In the case of a high-mass-loss-rate object such 
as CW Leo the line intensities change by < 5% and for a low- 
mass-loss-rate object such as R Lep by < 20%. Thus, no sig- 
nificant systematic trends are introduced through the choice of 
the SiO abundance description. Typically, we believe the abun- 
dances to be accurate to about + 50%, within the adopted cir- 
cumstellar model. 



3.5. Importance of vibrational excitation 

The excitation analysis includes the possibility of populating 
also the first vibrationally-excited state of SiO through the ab- 
sorption of 8 fim photons, primarily, from dust emission. This 
can significantly change the derived line intensities of transi- 
tions within the ground vibrational state as illustrated in Table|5] 
for CW Leo (and R Lep) for selected transitions. The calcu- 
lations were performed for a 15 m telescope. Neglecting to 
take account of the dust leads to about 50% lower integrated 
(over the line) intensities, which translates into at least a dou- 
bled fractional abundance of circumstellar SiO to produce the 
same flux. Interestingly, in this example approximating the dust 
emission using a single low-temperature (r = 510K) black- 
body containing the total luminosity of the source produces 
only about 10% higher intensities compared with the full model 
taking the dust emission/absorption fully into account. This is 



Table 5. Importance of IR pumping in the excitation of SiO. 





/(2 ^ 1) 


/(3 -. 2) 


/(5 ^ 4) 


/(8 ^ 7) 


Model with 


[Kkms^'l 


[Kkms^'l 


[Kkms^'l 


[Kkms '] 


CWLeo 










no (star + dust) 


11.2 


22.8 


41.1 


74.3 


star 


13.1 


28.0 


46.8 


81.0 


star + dust 


16.5 


47.3 


98.9 


159.2 


510KBB 


20.4 


51.5 


100.3 


172.6 


RLep 










no (star + dust) 


0.27 


0.60 


1.1 


2.0 


star 


0.34 


0.88 


1.7 


2.8 


star + dust 


0.34 


1.08 


2.6 


4.5 



not always the case and great caution should be taken when 
such approximations are made. It should be noted that CW Leo 
is a high-mass-loss-rate object and that the importance of the 
circumstellar dust radiation in relation to that of the stellar ra- 
diation field is less for a lower-mass-loss-rate object such as 
RLep(Table|5). 

4. Discussion 

4.1. Comparison with chemical models 

LTE stellar atmosphere models predict that the SiO fractional 
abundance in carbon stars is relatively low, typically ~ 5x 10"^ 
(see reviews by Glassgold 1999: Millar 2003. and references 
therein) about three orders of magnitudes lower than in M-type 
AGB stars. The relatively high SiO fractional abundances, de- 
rived in Sect. 13 of Ixl0"''-5xl0"^ can generally not be ex- 
plained by LTE chemistry. 

Departure from LTE could be caused by the variable na- 
ture of AGB stars that induces shocks propagating through 
the photosphere thereby affecting its chemistry. Models 
of shocked ca rbon-rich stellar atmospheres performed by 
IWiUacv & Cherchneff (1998) and Helling & Winters (20011) 
indicate that the SiO fractional abundance can be significantly 
increased by the passage of periodic shocks. There is also 
a strong dependence on the shock strength and possibly this 
mechanism can explain the observed fractional abundances, 
and their large spread, of SiO in carbon stars. 

From the results reported in Table|3it is clear that the frac- 
tional abundance of SiO varies substantially between the car- 
bon stars in the sample, from as low as 1 x 10"^ up to 5 x 10"^. 
As illustrated in Fig.0(solid circles) there is a clear trend that 
SiO becomes less abundant as the density, M/vg, in the wind 
increases. Even more intriguing, when compared to the distri- 
bution o f SiO fractional abundances for M-type AGB stars de- 
rived by |Gonz|lezDelgado.eial](|2003|) (Fig.0 filled squares), 
there appears to be no way of distinguishing a C-rich chemistry 
from that of an O-rich based on an estimate of the circumstellar 
SiO abundance alone. Also, the puzzling division of high-mass- 
loss-rate stars, into distinct low- and high-SiO fractional abun- 
dance groups, found bv Gonzalez Delgado et al. (2003) for M- 
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Fig. 6. Best-fit models (solid lines; parameters given in Table|5} for R Lep and R For overlayed on observed spectra (histograms). 
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Fig. 7. SiO fractional abundance (/o) obtained from an excita- 
tion analysis, as a function of a density measure (M/ve), for car- 
bon stars (filled circles) and M-type (O-rich) AGB stars (filled 
squares). The horizontal lines mark the abundances predicted 
from equilibrium chemistries. The dashed line shows the ex- 
pected f{oo) (scaled to 5x10"^, roughly the expected fractional 
abundance from stellar equilibrium chemistry when C/0<1) 
for a model including adsorpt ion of SiO onto dust grains (see 
lOonzalez Delgado et alJ2003[ for details). 



type stars is reenforced by the excitation analysis of the carbon 
stars. 

The non-detection of SiO maser emission (e.g., from 
the v=l, / = 2— » 1 transition) towards carbon stars (e.g., 
iLepine et al. 1978) is interesting in light of the high circumstel- 
lar SiO abundances found. Since SiO maser emission detected 



from M-type AGB stars emanates very close to the photosphere 
and well within the acceleration region (e.g., ICotton et alJ 
|2D04) this could suggest that in carbon stars the SiO molecules 
form somewhat further out in the wind where pumping by IR 
photons, and hence excitation of vibrationally excited states, is 
less effective. 



4.2. Condensation onto dust grains 

The importance of dust grains in regulating the circumstel- 
lar fractional abundance of SiO is supported by the observed 
trend that the SiO fractional abundance decreases with increas- 
ing density (M/Vg) of the wind (Fig.Q. This trend is expected 
from a classical treatment of the condensation process, as de- 
scribed in Gonzalez Delgado et al. (2003), where the conden- 
sation efficiency strongly depends on the dust mass-loss rate. In 
Fig.|3a depletion curve based on the simple model presented in 
[Gonzalez Delgado et al. (2003) is shown. Since the theoretical 
condensation results are very sensitive to the adopted parame- 
ters a large scatter around this curve is natural. 

In Fig.0the derived abundances are plotted against a den- 
sity measure where the terminal expansion velocity of the wind 
has been used. Since much of the condensation may occur in 
the acceleration region a lower velocity (and hence a higher 
density) is perhaps more suitable. However, if the character- 
istic velocity at which dust condensation is close to complete 
scales with the terminal expansion velocity (independently of 
the mass-loss rate), the results would only be shifted along 
the density axis in Fig. If, on the other hand, this veloc- 
ity is independent of the mass-loss rate plotting the abundance 
against the mass-loss rate is to be preferred. However, this only 
marginally changes the appearance of Fig.0(see also Fig. 5 in 
Gonzalez Delgado et al. 2003 ). 

The inclusion of adsorption of molecules onto dust grains 
is a shortcoming of current generation chemical models and 
will be required in a full model describing the chemistry of 
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AGB stars. In addition, grain surfaces could act as catalysts for 
chemical reactions. Although LTE chemistry can explain the 
SiO fractional abundances for a few of the carbon stars with 
the highest mass-loss rates we find it more likely that these 
abundances are the result of SiO depletion in the wind and that 
the photospheric abundance is significantly higher. High spatial 
resolution interferometric observations for a sample of sources 
could help to solidify this claim. 

4.3. The HCN/SiO line intensity ratio 

The HCN/SiO line intensity ratio has been shown to be a use- 
ful tool to discriminate between O-rich and C-rich envelopes 
(Buia rrabaletal. 1994; Olofsson et al. 1998; Bieging et al. 
I2OO0I) . The explanation for this has been that in O-rich en- 
velopes the abundance of HCN is very low and that of SiO 
high, naturally producing a low HCN/SiO line-intensity ratio. 
In a C-rich envelope HCN is abundant whereas SiO is not, pro- 
ducing a high HCN/SiO line-intensity ratio. Caution has to be 
taken here since optical depth effects are present, in particular 
for high-mass-loss-rate objects. 

From the present analysis it is suggested that the SiO frac- 
tional abundance in carbon stars closely mimic that in M-type 
(O-rich) stars. Therefore, we conclude that the observed dif- 
ference in the HCN/SiO line intensity ratio of about an order 
of magnitude is mainly due to a significantly lower fractional 
abundance of HCN in M-type AGB stars compared with the 
carbon stars. This claim has to be verified by performing a 
similar excitation analy s is, as presented for SiO here and in 
lOonzalez Delgado et alJ ( 120031) . also for a statistically signifi- 
cant sample of sources detected in HCN line emission. 



5. Conclusions 

New multi-transition millimetre SiO line observations of a 
sample of carbon stars are presented. With this addition, a to- 
tal of 19 carbon stars have been detected in SiO line emission. 
From a detailed excitation analysis, based on a reliable physical 
model of the sources, we reach the following conclusions: 

- In order to derive reliable fractional abundances of SiO, ex- 
citation through IR ro-vibrational transitions, needs to be 
taken into account. This will in most cases mean that the 
effect of dust grains need to be included in the excitation 
analysis. The present work contains a proper treatment of 
this. 

- The fractional abundance of SiO in carbon stars can be sev- 
eral orders of magnitude higher than predicted by thermal 
equilibrium chemistry. 

- A possible explanation for the high SiO fractional abun- 
dances found is a shock-induced chemistry. However, the 
influence of dust grains, both as a source for depletion as 
well as production of SiO, needs to be further investigated. 

- As observed for M-type AGB stars, a clear trend that the 
SiO fractional abundance in carbon stars decreases as the 
density (M/vg) of the wind increases is found, indicative of 
adsorption of SiO onto dust grains. 



However, these conclusions still rest on somewhat loose 
ground, e.g., the constraints on the SiO abundance distribu- 
tion are poor, and the relative importance of freeze-out onto 
dust grains, photodissociation, and circumstellar chemistry is 
still uncertain. More observations are required, in particular in- 
terferometric observations and high-/ single dish observations, 
to more firmly establish these claims. Also, a larger number 
of low-mass-loss-rate carbon stars must be observed with high 
sensitivity. 
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